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Abstract — Wireless networks can be self-sustaining by har- 
vesting energy from radio-frequency (RF) signals. Building on 
classic cognitive radio networks, we propose a novel method for 
network coexisting where mobiles from a secondary network, 
called secondary transmitters (STs), either harvest energy from 
transmissions by nearby transmitters from a primary network, 
called primary transmitters (PTs), or transmit information if 
PTs are sufficiently far away; STs store harvested energy in 
rechargeable batteries with finite capacity and use all available 
energy for subsequent transmission when batteries are fully 
charged. In this model, each PT is centered at a guard zone 
and a harvesting zone that are disks with given radiuses; a ST 
harvests energy if it lies in some harvesting zone, transmits fixed- 
power signals if it is outside all guard zones or else idles. Based 
on this model, the spatial throughput of the secondary network is 
maximized using a stochastic-geometry model where PTs and STs 
are modeled as independent homogeneous Poisson point processes 
(HPPPs), under the outage constraints for coexisting networks 
and obtained in a simple closed-form. It is observed from the 
result that the maximum secondary throughput decreases linearly 
with the growing PT density, and the optimal ST density is 
inversely proportional to the derived transmission probability 
for STs. 

I. Introduction 

Powering mobile devices by harvesting ambient energy such 
as solar energy, kinetic activities and electromagnetic radiation 
makes mobile networks not only environmentally friendly but 
also self-sustaining. In this paper, we propose a novel method 
for network coexisting where mobiles in a cognitive-radio 
(secondary) network harvest radio-frequency (RF) energy from 
transmissions in a primary network besides opportunistically 
accessing the spectrum licensed to the primary network. The 
throughput of the secondary network is studied based on a 
stochastic-geometry model where the primary and secondary 
transmitters (PTs and STs) are modeled as independent homo- 
geneous Poisson point processes (HPPPs). In this model, each 
PT is protected from the interference from STs by a guard 
zone and delivers significant RF energy to STs located in a 
harvesting zone. This model is applied to maximize the spatial 
throughput of the secondary network over the secondary 
transmission power and node density under outage-probability 
constraints for the primary and secondary networks. 

Cognitive radios enable efficient spectrum usage by allow- 
ing a secondary network to share the spectrum licensed to a 
primary network [1]. Besides active development of algorithms 
for opportunistic transmissions by secondary users (see e.g., 
[2], [3]), research has been carried out on characterizing 



the throughput of coexisting networks based on stochastic- 
geometry network models. The trade-off of the capacities 
between two coexisting networks [4], [5] and among multiple 
networks [6] has been studied. Moreover, outage probability 
of a Poisson-distributed cognitive radio network considering 
guard zones has been analyzed in [7]. Our work also considers 
Poisson-distributed cognitive radio networks but differs from 
the prior studies in that the secondary users power their 
transmitters by harvesting energy from primary transmissions 
instead of using traditional power supplies such as batteries. 

Recently, wireless communication powered by energy har- 
vesting has emerged to be a new and active research area. 
Given energy harvesting, the existing transmission algorithms 
need to be redesigned for different wireless systems [8], 
[9] to account for the randomness and causality of energy 
arrivals. For broadcast channels, multi-antenna transmission 
for simultaneous wireless information and energy transfer has 
been studied in [10]. Also, the throughput of a MANET 
powered by energy harvesting has been investigated in [11] 
where the network spatial throughput is maximized by op- 
timizing transmission power under an outage constraint. In 
contrast, our work studies a different type of network, namely 
cognitive radio network powered by opportunistic RF-based 
energy harvesting [12]. 

In this paper, it is assumed that time is slotted, and both 
the channels and HPPPs are quasi-static. The harvesting 
zone and secondary transmission power are constrained to 
be sufficiently small such that a ST lying in a harvesting 
zone is fully charged within a time slot. The main contri- 
butions of this work are summarized as follows. Based on 
the aforementioned network model, we derive the transmission 
probability of STs using Markov-chain theory. Moreover, the 
maximum network throughput of the secondary network under 
the outage constraints is derived by a Poisson approximation of 
the transmitting STs. It is shown that the maximum secondary 
throughput is linearly decreasing with the increasing PT den- 
sity, and the optimal ST density is inversely proportional to 
the transmission probability of STs. 

II. Model and Metric 

A. Network Model 

As shown in Fig. 1, the PTs and STs are distributed as 
independent HPPPs denoted by $ p = {A"} and $ s = {Y} 
with density A p and A s , respectively, where X, Y e R 2 
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Fig. 1: A stochastic -geometry model of the cognitive radio 
network where PTs and STs are distributed as independent 
HPPPs. Each PT/ST has its intended receiver at an unit 
distance (not shown in the figure). A ST inside the harvesting 
zone harvests energy from the transmission by a nearby PT. 
To protect primary transmissions, all STs inside guard zones 
are prohibited from transmission. 

denote the coordinates of the PTs and STs, respectively. Each 
transmitter transmits data with fixed power to an intended 
receiver at an unit distance. Time is partitioned into slots with 
unit duration. The point processes are fixed within each slot 
and vary independently over different slots. STs are prevented 
from transmission when they lie in the guard zones to protect 
primary transmissions. Specifically, consider a disk with radius 
r g centered at each X and let 



g= |J b(X,r g ) 



(1) 



xe<s> v 



where b(x, r) C H 2 represents a disk of radius r centered at 
x and hence b(X, r g ) is the guard zone with radius r g for 
protecting X E <£> p . A ST 7 e $ s cannot transmit when 
Y E Q. Note that intuitively a guard zone should be centered at 
a primary receiver but is defined as above to simplify analysis 
[2]. Let R denote the distance between a typical ST, denoted 
as Y', and its nearest PT. The probability p g that Y' lies in a 
guard zone is given as [13] 



Pg = Pr(y' e Q) 

= Pr(i? < r g ) 
= 1 - e -* X r r *. 



(2) 
(3) 
(4) 



B. Energy-Harvesting Model 

Each ST Y E Q s harvests energy from transmission by 
the nearest PT when it is located inside a harvesting zone. 
Let b(X, r/j) with < < r g represent the harvesting 
zone centered at X E & p such that a ST Y E b(X, ru) 
harvests energy from the PT X. Channels are assumed to 
have no fading for simplicity. Let P p and P s denote the fixed 
transmission power of PTs and STs, respectively. It is assumed 
that P p 3> P s , Given R < r^, the harvesting power of Y' is 
T]P p R~ a where < r\ < 1 denotes the harvesting efficiency 
and a > 2 is the path-loss exponent. Similarly as (4), the 



probability p^ that a typical ST lies inside a harvesting zone 
is 

(5) 



Ph = l-e-*** r ». 



We assume that each ST has finite battery capacity P s that 
is the same as its transmission power. Upon the battery being 
fully charged, a ST transmits in next slot if it is outside all 
guard zones. 

C. Performance Metric 

Let p t denote the transmission probability of a typical 
ST Y'. Assuming interference-limited networks, the received 
signal-to-interference-ratio (SIR) is required to exceed a target 
SIR for reliable transmission. Let 9 P and S be the target SIRs 
for the typical primary and secondary receivers, respectively. 
The outage probability is then defined as P^ = Pr(SIR (p) < 
6 P ) for the primary network and P ( $ = Pr(SIR (s) < 0„) 
for the secondary network. Outage-probability constraints are 
applied such that < e p and P^( < e s with < e p , e s < 
1. Given fixed PT density A p and transmission power P p , 
the performance metric is the spatial network-throughput C s 
(bit/s/Hz/unit-area) of the secondary network given by 



C s = p f A s log 2 (l + 9 S ) 



(6) 



under the outage constraints. 



III. Transmission Probability of Secondary 
Transmitter 

In this section, the transmission probability of a typical ST 
is analyzed. Note that the minimum power harvested by a 
ST is r)P p r^ a . Therefore, the battery of an energy-harvesting 
ST is fully charged within a slot if < P s < T]P p r~^ a and 
this scenario is assumed in this paper. It then follows that 
the battery power level can only be or P s . Consider the 
finite-state Markov chain with the state space {0, 1} and let 
the state and 1 denote the battery level of power and 
P s , respectively. Furthermore, let P represent the transition- 
probability matrix that can be obtained as 



1 ~Ph 

l-Po 



Ph 



(7) 



with pg and ph given in the preceding section. The steady- 
state probability is then given by the vector n = [no 7Ti] using 
Markov-chain theory as shown in the proof of the following 
proposition. Note that 7Ti represents the probability that the 
battery is fully charged. 

Proposition 1. For < P s < r\P p r~^ a , the transmission 
probability of a typical ST is 



Pt 



(8) 



PhQ- - Pg) 
Ph + l-Pg' 

Proof: TV is the left eigenvector of P corresponding to 
the unit eigenvalue such that 



ttP 



(9) 



From this equation, the distribution of the battery level at a 
typical ST is obtained as 



1 



7T = 



Pg 



7T1 



Ph 



Ph 



1 - 



(10) 



Consequently, the desired result in (8) is obtained by multi- 
plying 7Ti with 1 — p g . ■ 
It can be observed from the above result that p t is inde- 
pendent with P s , provided that < P s < rjP p r^ a , since the 
battery of an energy-harvesting ST is fully charged over one 
slot if it is inside a harvesting zone. 

IV. Network Throughput Maximization 

In this section, the network throughput of the secondary 
network defined in (6) is investigated under the outage con- 
straints. 

A. Outage Probability 

There are four types of interference in the network model. 
Let I pp , I ps , I ss and I sp denote the aggregate interference 
power from all PTs to a typical primary receiver/secondary 
receiver, and from all transmitting STs to a typical secondary 
receiver/primary receiver, respectively. Then the outage prob- 
ability for the typical primary receiver located at the origin 
can be written as 



,(p) 



where I pp = 



Pr 



Inn ~\~ Is 



< 



and \X\ € 



(11) 



is the distance 



between X and the origin. It is important to note that I sp 
is a summation over the process of transmitting STs that 
is not necessarily a PIPPP. As a result, the analysis of the 
distribution of I sp is intractable. To overcome this difficulty, 
we approximate the transmitting ST process as a HPPP with 
density p t \ s and assume that it is independent with $ p . This 
assumption is validated by simulation in the sequel. It is worth 
mentioning that a similar approximation is also used in [7]. 
Based on this approximation, we obtain the following lemma. 

Lemma 1. By approximating the process of transmitting STs 
as a HPPP with density Pt\ s , the outage probability of the 
primary network is given as 



(12) 



where A(y) and T £ M 2 denote a HPPP with density y and 
the coordinate of a node, respectively, and 




Tp = dp Pt^s 



P, 



An 



(13) 



Proof: Let $ t C $ 
transmitting STs. Then L 



sp 



Er e $ t \t(x«,r s ) p s\ Y \ a where X* is the PT corresponding 
to the typical primary receiver at the origin. Let /' = 



follows that 



P S \T\ a and we approximate I sp as I' Then, it 



J 4- V 



Pr 



= Pr 



]T \T\~ a + ]T \T\~ a >l 

\TeA(\[ p> ) TeA(A< p) ) 

/ 



E 



\T\~ a > 1 , 



(14) 



(15) 



(16) 



\TeA(\[ p) +\i p) ) 



where = Q 2 J a \p, = {^p^ Pt^s, and hence 

A^ p) + \ { 2 p) = t p . Note that (15) is from the fact that the point 
process aA(A) with a > is a HPPP with density A/a 2 which 
can be proved by the Mapping Theorem [11], and (16) is from 
the Superposition Theorem [13]. ■ 

It is worth mentioning that P^ depends only on t p and is 
an increasing function of t p . 

For the secondary network, it is important to note that the 
outage probability for the typical secondary receiver located 
at the origin is implicitly conditioned on that the typical 
ST corresponding to this receiver is active (transmitting). 
Consequently, the typical ST must be lying outside guard 
zones and hence there are no PTs inside the disk of radius 
r g centered at the typical ST. Let the event of this condition 
be denoted by £ = {$ p n b(Y*,r g ) = 0}, where Y* denote 
the transmitting ST corresponding to the receiver at the origin. 
Then P^) can be written as 



Pr 



Ps 



where I ps — I pp and 



<0 S \£ 



P,\Y\ 



(17) 



From the law 



of total probability, 

Pr( 7 * 



P W _ 

^out — 



+Is 



<^)- Pr ( I ^fc <0s\£)Pr(£) 



(18) 

Note that £ = {$ p n b(Y*,r g ) ^ 0}. Similar to Lemma 1, 
I' sp is a good approximation of / ss , and hence we obtain the 
following lemma. 

Lemma 2. By approximating the process of transmitting STs 
as a HPPP with density PtX s , the outage probability of the 
secondary network can be approximated as 



P W 

■'out 



Pr 



TgA(t s ) 



\T\ 



> 1 



Pg 



1 



Pg 



denote the point process of where 
can be written as I sp = 



t s = s a p t \ s + X p 



(19) 



(20) 



Proof: Assuming that I' is a good approximation of I ss , 
it then follows that 

^(r^r < *.) - ^ (r4r^ < 9s ) (21) 



Pr 



\T\~ a > 1 , (22) 



where A 



yTeA(A^ s '+A!j 

s 2/ > t A s , A 2 s) 



-2/a 



A , and thus 

(s) (s) 

X\ ' +A 2 = r s- Note that (22) follows the same arguments as 
for (15) and (16). Next, with the aforementioned assumption 
that P p ^> P s , it is likely that interference from even only 
one PT inside b(Y*,r g ) is sufficient to cause outage at the 
typical secondary receiver at the origin. Consequently, we have 



Pr 



Op 



the fact that Pr(£) 
result in (19). 

Note that similarly as P out 
of To, 



1. Substituting this result, (22) and 
1 — p g into (18) yields the 



e - nX » r 9 



(p) p( s 



P out is an increasing function 



B. Network Throughput Maximization 

The throughput of the secondary network, defined in (6), is 
maximized over P s and A s , for a given pair of P p and X p and 
under primary and secondary outage-probability constraints. 
The optimization problem can be written as 



(PI): max p t A s log 2 (l 



+ S ) 



s.t. 



pO) < , 

p(«) < 



(23) 

(24) 
(25) 



As stated in Section III, we assume < P s < rjPpr^ 01 in this 
paper and thus the transmission probability p t is a constant for 
a given pair of r g and r^. Therefore, the maximum network 
throughput is achieved by simply maximizing A s under the 
outage constraints. To solve Problem (PI), we use the nominal 
node densities [11], denoted as /i p and p, s , for HPPPs A(/z p ) 
and A(/x s ), respectively, such that 




T\~ a > 1 = (l-p g )e s +p g . 



(26) 



(27) 



It should be taken into account that jj, p and p, s are monotone- 
increasing functions of e p and e s , respectively, which in gen- 
eral have no closed-form expressions and can only be found 
by simulation. With these parameters, the solution of Problem 
(PI) can be found as given in the following proposition. 

Proposition 2. The maximum network throughput of the 
secondary network is given by: 





Optimal point 


A, = f,(P, 




Admi^sitjle set 


\, = h{P, 








y p V J p V f^p j 

Fig. 2: The admissible set of (P S ,X S ) shown as the shaded 
region for Problem (PI). 



(28) 



CI = (97% fi, - V-°r 2 h \ p ) log 2 (l + 67), 
where the optimal ST transmit power is 

P: = »W> 
and the optimal ST density is 



(29) 



(30) 




lis 



where the optimal ST transmit power is 



X p log 2 (l + s ), (31) 



P* = 



fts 



p \ ftp 

and the optimal ST density is 



p 



Pt 



Js "lis 



I'- 



lh' 



(32) 



(33) 



Proof: Since P~$ and P ^t 316 increasing functions of t p 
and t s , respectively, Problem (PI) is equivalent to maximizing 
A s subject to T p < Hp, t s < fi s and < P s < ?yP p r^ Q . 
In addition, from (13) and (20), t p < /j, p and r s < fi s are 
equivalent to X s < fi(P s ) and X s < /2(Ps), respectively, 
where 



h{p s ) = - (e p 

Pt v 



f2(P.) 




(34) 
(35) 



As plotted in Fig. 2, fi(P s ) decreases whereas /a(-Ps) in- 
creases with growing P s . The intersection point of the two 
curves X s — fi(P s ) and X s — /2(-P) can be found 
by solving fi(P s ) = /a (-Pa), which yields (P S ,X S ) = 

fJ> s (fJ>p — 6p A p ) 



^ ^ p 



pi 



The shaded region in Fig. 2 

shows the admissible set of (P s , A s ) that satisfies A s < f%(P s ), 
A s < fiiPs) and < P s < rjP p r^ a . It can be observed 
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Fig. 3: Comparison between the approximated and simulated 
outage probabilities with \ p = 0.01 and X s = 0.1. 



from Fig. 2(a) that if t]P p r h 



< 



9* (BsX 



P.„, \ s attains 



its maximum when P* = T)P p r h a . Thus, the results in (28) 
and (30) are obtained from A* = / 2 (P*). Otherwise, if 

vP P rh a > t (fA 2 Pp> it is observed from Fig. 2(b) that 
the optimal value of A s is the intersection of X s — fi(P s ) and 
= f2{P s ), as obtained above. ■ 

Note that if r/r^" < |*- \£~J > the secondary trans- 
mission power is small enough such that satisfying the sec- 
ondary outage constraint guarantees satisfying the primary 
outage constraint. Also note that C* for the case of r)r^ a > 

k fe) 2 is alwa y § lar s er than that of ^h a < e t v fe) 2 ■ 

Two remarks on the trade-offs between the two coexisting 
networks are in order. 

1) The maximum secondary network throughput C* de- 
creases linearly with growing PT density A p . Also, the 
optimal ST density A* is inversely proportional to the 
secondary transmission probability pt. 

2) If A p goes to zero, C* converges to its largest value; 
however, since p t converges to zero as A p — > 0, 
A* diverges to infinity. It means that although sparse 
PT density results in larger secondary throughput, an 
extremely large number of STs should be deployed to 
achieve the maximum value. 

V. Numerical Results 

Numerical results are provided in this section. The path- 
loss exponent is set as a = 4, the power ratio between the two 
networks is If = 0.1, the harvesting efficiency is ri = 0.1, and 
the radius of guard zone, harvesting zone is r g = % ry l = \. 

Fig. 3 compares the outage probabilities obtained by simu- 
lation and those based on approximations in (12) and (19). It 
is observed that our approximations are quite accurate. 

The relation between A* and X p for the case of ?77\ Q > 

I s - (j^j is plotted in Fig. 4(a) with three different values 
of e p and e s = 0.4. As mentioned above, A p —> results 
in deploying an infinitely large number of STs to achieve 
the maximum secondary network throughput. Moreover, it is 




(a) Optimal ST density (b) Maximum secondary throughput 

Fig. 4: Optimal ST density A* and maximum secondary net- 
work throughput C* versus PT density A p with 8 p = 8 S = 5. 

observed that larger e p allows deploying more STs as it permits 
more interference from the secondary network. 

The maximum secondary network throughput C* versus A p 
is plotted in Fig. 4(b). The parameters are set as e p = 0.35 

and e s = 0.4 for the case of r)r^ a < |j (fr) 2 ' £ p = 

0.2 and e s = 0.4 for the case of r)r^ a > ^ (^j 2 . It is 
observed that C* decreases linearly with increasing A p . 

VI. Conclusion 

In this paper, a novel scheme is proposed to enable the 
secondary network to harvest energy and reuse the spectrum 
in the primary network. The secondary network throughput is 
maximized and derived in a closed form based on a stochastic- 
geometry network model under outage constraints. 
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